APS/123-QED 



Multimode analysis of the light emitted from a pulsed optical parametric oscillator 

Anne E. B. Nielsen and Klaus M0lmer 
Lundbeck Foundation Theoretical Center for Quantum System Research, 
Department of Physics and Astronomy, University of Aarhus, DK-8000 Arhus C, Denmark 

(Dated: February 1, 2008) 

We present a multimode treatment of the optical parametric oscillator, which is valid for both 
pulsed and continuous-wave pump fields. The two-time correlation functions of the output field are 
derived, and we apply the theory to analyze a scheme for heralded production of non-classical field 
states that may be subsequently stored in an atomic quantum memory. 
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I. INTRODUCTION 

The process of parametric down conversion, in which 
pump photons with frequency u> are sent through a non- 
linear crystal and converted into pairs of signal and idler 
photons with frequencies w\ and u>2 = w — u>i, is a well- 
known and widely used phenomenon in quantum optics. 
Since the conversion efficiency is rather small, in the opti- 
cal parametric oscillator (OPO) the process is enhanced 
by placing the crystal inside an optical cavity. Due to 
its importance as a source of squeezed light when oper- 
ated below threshold, the OPO has been subject of much 
investigation, in particular for the case of a time indepen- 
dent continuous- wave pump field [l|, [H, Q . The case of 
a periodically modulated pump field has been studied in 
[J]. In the standard treatment of the OPO, the calcula- 
tions are performed in frequency space, the pump is taken 
as a time independent monochromatic beam, and it is of- 
ten assumed that the cavity allows only one or two modes 
with well-defined frequencies, depending on whether de- 
generate or nondegenerate operation is considered. The 
field leaking out of the cavity is then determined from 
input-output formalism. 

An OPO may also be pumped with a pulsed light field. 
This case is often discussed as if the OPO generates a 
single mode pulse of squeezed light, but this is not quite 
the case. The output consists of a number of independent 
modes squeezed by different amounts as discussed for sin- 
gle pass down conversion in Refs. 0, In the present 
paper we use a completely different approach to charac- 
terize the output from a pulsed OPO. The treatment is 
a generalization of the discussion of the continuous- wave 
OPO in Ref. [7[, and in contrast to Refs. [ll, S S S la| all 
calculations are performed in the time domain. 

The time domain treatment has several advantages. 
Firstly, it is not necessary to assume the existence of 
a single mode or a couple of independent modes in the 
cavity since all effects of cavity resonances appear natu- 
rally in the analysis from the specification of the length 
of the cavity via the cavity round trip time. Further- 
more, since we do not assume orthogonal cavity modes, 
our treatment is valid for all values of the transmission of 
the OPO output coupling mirror and not just for small 
transmissions. In particular, we may increase the mirror 



transmission to unity in our formulas and obtain the re- 
sults for single pass down conversion. Secondly, we do not 
assume any particular temporal shape of the pump field, 
and the analysis is thus suitable to investigate the transi- 
tion between the few mode situation for short pulses and 
the highly multimode regime for continuous-wave fields. 
In the limit of a time independent pump field our re- 
sults reduce to those of Ref. [l| , provided we invoke extra 
approximations which guarantee the existence of only a 
single cavity mode. Thirdly, the time domain analysis is 
convenient if the light emitted from the OPO is detected 
continuously in time and the back action of the measure- 
ments on the system has to be taken into account, see 
0. 

Although the OPO output can be both squeezed and 
entangled, it belongs to the family of so-called Gaussian 
states, and for a number of applications in quantum com- 
munication and computing, non-Gaussian states are nec- 
essary. The OPO output can be transformed into non- 
Gaussian states by measurements as shown schematically 
in Fig. [TJ a small fraction of the light is extracted, pos- 
sibly frequency filtered, and detected by an avalanche 
photo diode (APD). As, e.g., a squeezed vacuum state is 
a superposition of even photon number states, a photo 
detection in the APD heralds the generation of a state 
with only odd photon number contributions such as sin- 
gle photon states and odd Schrodinger kitten states in the 
remaining beam as verified experimentally in 0, [Io| . 
Ref. [ll[ provides a theoretical determination of the mode 
of the output state with the largest single photon fidelity 
for the case of a continuous- wave nondegenerate OPO. 
We note that in the experimental verification by homo- 
dyne detection, one may use a constant amplitude local 
oscillator, measure the homodyne detector signal as a 
function of time, and then at a later stage multiply the 
recorded data with the mode function, selected according 
to the instant of the photo detection by the APD [9( . 

For several applications of the heralded non-Gaussian 
field states it is, however, necessary to mode match the 
state in real time. This is the case if one wants to interfere 
different pulses and if one wants to store the field state in 
an atomic medium and, e.g., perform quantum gates by 
optical non-linearities in the medium [12j . Whether mak- 
ing use of electromagnetically induced transparency [l3[ , 
a Raman transfer [141 ] . or Faraday polarization rotation 
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FIG. 1: Experimental setup to generate nonclassical states of 
light. BS denotes a beam splitter with a small reflectance, 
and the curve to the right is an example of a possible mode 
function that may contain a non-Gaussian state conditioned 
on a detector click in the APD. 

[l6| , storage of light pulses in atomic media involves 
a strong pump pulse that defines the temporal mode of 
the quantized field to be mapped onto the atomic col- 
lective degree of freedom [l7|. This implies that storage 
of a field state heralded by an APD detection event re- 
quires a storage control field, which has to be produced 
at the appropriate time. Since the mode function oc- 
cupied by the non-Gaussian state may extend to times 
both before and after the trigger event [ll|, it may thus 
be necessary to switch on the control field before the trig- 
ger event or, more realistically, to delay the arrival of the 
signal to the storage medium. For demonstration exper- 
iments, it would be attractive to apply the storage pro- 
tocol irrespective of the APD output and subsequently 
post select the instances where APD detection actually 
occurred in a suitable time window around the instant 
for which the stored mode is optimal. In this paper we 
do not model the storage process itself, but we identify 
the mode functions that maximize the non-Gaussian fea- 
tures of the pulses to be stored, conditioned on realistic 
detection events. In particular, we wish to identify mode 
functions that are optimal independent of the precise de- 
tection time within a predefined time interval, which can 
not be chosen arbitrarily small if one wants a satisfac- 
tory APD detection probability. We shall for concrete- 
ness take the optimal mode to be the one with the most 
negative value of the Wigner function at the origin. 

The paper is structured as follows. In Sec. II we de- 
velop a theoretical description of the OPO and derive the 
two-time correlation functions of the OPO output field. 
In Sec. Ill we show that a treatment similar to the one 
applied to the OPO can be used to express the output 
field from the filter in Fig.[T]in terms of the input field. In 
Sec. IV we determine the mode of the conditioned state 
with the most negative value of the Wigner function at 
the origin, we compute the probability to obtain a trigger 
detection event within a specified time interval, and we 
compare the performance of continuous- wave and pulsed 
operation of the OPO for preparation of pulses that can 
be stored with predefined control pulses in an atomic 
medium. Section V concludes the paper. 



II. CORRELATION FUNCTIONS FOR THE 
OPTICAL PARAMETRIC OSCILLATOR 
OUTPUT FIELD 

A theoretical model of the OPO is based on the setup 
illustrated in Fig. [5] We have chosen a four-sided ring 
cavity for mathematical convenience, but this is not es- 
sential to the analysis. The field annihilation operators 
of the input and the output fields are denoted a(t) and 
b(t) 1 respectively, while Ci(t), £ = 1,2, represents the field 
at different positions inside the cavity as shown, and v(t) 
is the annihilation operator of a field in the vacuum state. 
The beam splitter BSi couples the input and output fields 
to the intra cavity field, and the fictitious beam splitter 
BS2 models losses in the system. We model the entire 
loss as if it takes place between the crystal in the upper 
part of the figure and the output coupling mirror, which 
represents a worst case situation. 

The crystal is pumped by a classical pump field 
f(t) = |/(f)|e^^'\ which is assumed to pass unhindered 
through the cavity mirrors. The parametric process in 
a non-linear crystal is described in [g, |6(, and leads, in 
the time domain, to the following mapping of the field 
incident on the crystal to the field leaving the crystal 

c(t) -> cosh(2 X |/(i)|)c(<) - ie^'W sinh(2 X |/(i)|) C t(i). 

(1) 

In (fljl y. which is taken to be real, is proportional to 
the second order susceptibility of the crystal and to the 
length of the crystal. We have assumed that the band- 
width of the down conversion is infinite and neglected 
depletion of the pump field and differences in the phase 
matching conditions for down conversion to different fre- 
quencies. This is likely to be a good approximation in the 
experiment described in Ref. [9|], since the crystal inside 
the OPO is only few millimeters long, and thousands of 
cavity modes are populated. Also, modes far from the 
center of the spectrum are irrelevant because they are 
filtered out by other components in experiments (see the 
next section). 

In order to express b(t) in terms of a(t) and v(t), we 
use the beam splitter relations 

b(t) = hd 2 (t) + ina(t) (2) 

CiO) = i ri c 2 (t) + ha(t) (3) 
and express c 2 (t) in terms of c\ (t — r) 

c 2 {t) = — it<x ( cosh (z(t — t)) c±(t — t) 

- ze i0(t ~ r) sinh (z(t - t)) c\(t - r)) + ir 2 v{t). (4) 

ti (rt) is the field transmission (reflection) coefficient of 
BSi, and r is the round trip time in the cavity. In ij^j). 
we have assumed that each cavity mirror induces a phase 
shift of 7r/2, and we have defined the following quantities 

z(t-r) = 2x|/(t-r + fr)| (5) 
</>(t-r) = 0/(*-r + fr), (6) 
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FIG. 2: Model of an optical parametric oscillator. Photon 
pairs are generated by parametric down conversion in the 
crystal, which is pumped by the classical field f(t). The beam 
splitter BSi couples the intra cavity field to the input (vac- 
uum) field a(t) and the output field b(t), and the beam splitter 
BS2 is included to model intra cavity losses. 



where £r is the traveling time from BSi to the crystal. 
Furthermore, for convenience, we have redefined v(t) ac- 
cording to v(t — C T ) ~^ w(t), where £t is the traveling 
time from BS 2 to BSi. Isolating b(t) from Eqs. ©, ©, 



and ((4} and assuming e 



#(*) 



T ) for all t, we find 



S(i) =ir 1 a{t)-it 2 1 Y J r i t 2 +1 

( cosh ^2 z ( f - fcr )^ «(* ~{n+ 1)t) 
- ie lc/ > {t) sinh (J2 z ( f ~ fcT ) « + (*- (» + l) r ) 

OO 

+ ihr 2 v{t) + itxr 2 ^(rii 2 )" +1 

n=Q 

( cosh z(t - fcr)^ u(f - (n + l)r) 
- ie^ ( *' sinh (j2 z (t ~ kr ) J V f (* - (« + !) r ) J • 



(7) 



vfc=l 



The requirement e 4 "^*' = e l ^ l ~ T \ which is satisfied for 
(j){t) = — 2irNt/r + (f>o, where N is an integer and <po is a 
constant, means that the successive squeezing operations 
add up in phase, i.e., it is the same quadrature that is 
squeezed at times t, t ± t, t ± 2r, .... 

For the special case of a time independent continuous- 
wave pump field z{t) = z, below threshold r^e 2 < 1, a 
Fourier transform of ([7|) leads to the following expression 
for the output field in frequency domain: 

6(w +^) = Gi{uj+LjQ)a{uJa+ijj)+gi(oj~oj a )a) {ujq — uj) 
+ G 2 (io + lo )v(lo + u) + ff2(w - uo)v'(wo - uj), (8) 



where loq — Nt:/t is half the frequency of the pump field, 



Gi(lo + lo q ) = in - i-^-G+{u + w ), (9) 
G 2 (w + w ) = «tir 2 +i G + (w + w ), (10) 



-e^!^ G _ (w - Wo ), 



g 2 (io~oo ) = g G_(«-wo), 



and 



1 - ri^ 2 e z+l ' 



± 



1 - rii 2 e- z+i 



(11) 
(12) 

(13) 



Assuming that the input field a(t) is in the vacuum state, 
the frequency correlation function 

(P(w)b(uj')) = t\t\(\ - r\t\) sinh 2 (z)<5(o; - J) I 
1 + r\t\ + Ar\t\ cosh 2 (z) + 2r 1 t 2 (r 1 t 2 cos(2lot) 

-2(l + r 2 t^)cosh(z)cos(wr))) (14) 



shows that resonances occur for lot = n2ir, n G Z, as 
expected. In particular, the degenerate frequency loq is 
only a resonance frequency if N is even. To compare Eq. 
([5]) to the results given in Ref. [TJ] for a cavity with a single 
resonance frequency at loq, we assume N to be even and 
lot <C 1 in Eq. ([8]) . In order to have well separated cavity 
modes we must also assume t\ <C 1 and r| <C 1, but to 
maintain a finite width of the modes as in Ref. [l[ we 
keep 71 = tf/r and j 2 = t|/t fixed. Since tji <C 1 and 
t"/ 2 <C 1 in this limit, the mean number of round trips in 
the cavity is large, and z must be correspondingly small, 
i.e., z must be proportional to r. For \z\ = |e|r we obtain 
Eq. 46 of Ref. fjj. 

Returning to the general case, we compute the two- 
time correlation functions of the output field for a(t) in 
the vacuum state from Eq. ([7]): 



(b\t)b(t'))=tltt(l-rltl) 



00 00 



'm+l 



J2 int 2 ) q+2m smh ij2 kT ) 

, g— 1 m— 



, k=l 



/q+m+l \ 

sinh I ^ z(<' - fcr) j 5(i - 1' + gr)- 



OC OO 



/g+m+1 N 

EE(^) 9+2m smh 

q=0 m =0 \ fe=l / 

/m+l \ 
sinh ^ z(t' - fcr) (5(i - t' - qr) 



(15) 
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and 



T5l) and ([TBI) reduce to 



(S(t)S(f)> = ie^^Cl - r?t|) 



oo oo 



'm+1 



E J2 ( r ^) q+2m sinh ( £ z ^ ~ fcr ) 

, g— 1 m— 



, fe=l 



cosh - kr) S(t -t' + qr)+ 

oo co /q+m+1 \ 

^^(rit2) ?+2ra cosh z (t- k r) 

q=0 m=0 \ fe=l / 

sinh [ ^ - kr) ] 5(f - i' - gr) ] 



(16) 



, fe=i 



These two expressions and (b (t)} — are sufficient to 
characterize the output state completely, because it is 
Gaussian. 

The correlations of an arbitrary single mode with mode 
function h{t) and annihilation operator 



b= / h*(t)b(t)dt 



(17) 



are easily calculated from Eqs. (fT5|) and (fT6|) . Wc 
shall only be concerned with the degenerate case be- 
low, and we thus choose h(t) = fr(£) e ^o/2-M7r/4 e -i Wo t = 
~ h (tyi<t>{t)l2+iKl± with ^ real The constant phase fac _ 

tor is included to obtain a real value of (6 2 ), which cor- 
responds to the case, where the axes of the squeezing 
ellipse in phase space lie along the quadrature axes. In 
the following we take hit) to be real and omit the factor 
ie i4>(t) i n Eq. (fiTj]) . 



1 ~rj4 c (2+\g\)z ! l-rft 



1 - r\t\e 2z 



2+ 2 e -(2+| 9 |) 2 



and 



e |g|« _ e -\q\z 



2+2 



1 - r 2 i 2 e - 2z 

8(t-t'-qr) (18) 



<6(i)S(t')> = ie^>^ £ (^) 

q— — oo 
= (2+|«|)2 _ 



1 - „(2+| |)z 1 _ r l*2 



1 - r 2 i 2 e 22 



1 - r\t\e- 2z 



-(2+|«|)2 



e l9l* + e -l«l* Wi-i'-gr). (19) 



Until now we have described the fields in terms of time 
dependent Heisenberg picture operators, but it is also 
useful to consider the OPO model from a Schrodinger 
picture point of view. In particular, the latter approach 
is suitable for a numerical treatment of the pulsed OPO. 
To this end we divide the light beams into small segments 
of (infinitesimal) duration At and treat each segment as 
a single mode, i.e., we define 



,.s _ / 1/y/Ai for ti - At/2 <t<U + At/2 
9l{ ' ~ I otherwise 



(20) 



di = / g*{t)d(t)dt = d{U)VAt 



The variances of x = ib+V) / \/2 andp = —i(b~tf)/V2 
are 1/2 + (fo+fo) + (b 2 ) and 1/2 + (6+fo) - (6 2 ), respectively, 
and the most efficiently squeezed mode is thus obtained 
by choosing h(t) as the eigenfunction with the smallest 
eigenvalue A of an integral equation with either the in- 
tegration kernel 5{t - t')/2 + (tt(t)b(t')) + (b(t)b(t')) or 
S(t - t')/2 + (tf(t)b(t')} - (b(t)b(t')), and A is the corre- 
sponding variance. Since the output field at time t is only 
correlated to the output field at times t + jit, n G Z , the 
mode function of the most efficiently squeezed mode is 
in general a spike function that is only nonzero at times 
t + tit, n = ...,—1,0,1,2,... , for some t, but if the vari- 
ations of the pump field take place on a time scale r or 
more slowly, displaced spike functions, and therefore also 
a smooth mode function with the same envelope, will be 
only slightly less squeezed. 

In the limit of a time independent pump field, Eqs. y = {x\,p\,x-i,P2 



and replace the continuous field annihilation operator 
d(t) (d(t) = d(t), b(t), ci(t), or C2(t)) by the discrete 
annihilation operators 



(21) 



localized at time ti. In practice, we will have to deal with 
a finite At, and neglecting the field variation within each 
At interval is equivalent to a cutoff in frequency, justi- 
fied by the experimentally relevant frequency regime. A 
thorough discussion of continuous and discrete operator 
descriptions of light beams in both time and frequency 
space may be found in Ref. [18j. Since the output field 
and the cavity field are Gaussian for a vacuum input 
field, the state of all the small light beam segments is 
efficiently represented by a Wigner function. In general, 
the Wigner function of an n-mode Gaussian state with 
zero mean values is on the form 



W(y) 



1 



■K n sJdet(V) 



(22) 



where y — (x\,p±, X2,P2, ■ ■ ■ , x n ,Pn) T is a column vec- 
tor of quadrature variables, and V is the covariance 
matrix of the n modes. V = (yy T ) + (yy T ) T , where 
,x n ,p n ) T , x l = (di + 3[)/yj2, 



pi = —i(di — dl)/\/2, and di is the field annihilation op- 
erator of mode i. Before the pump pulse reaches the 
crystal, all modes are in the vacuum state, and V is the 
identity matrix. We include all the r/At cavity modes 
and a sufficiently large number of a(t) and v(t) modes in 
(|2"2")l . As time passes by modes are squeezed, when they 
hit the crystal, input and cavity modes are transformed 
into cavity and output modes, when they hit BSi, and 
vacuum and cavity modes are transformed into cavity 
and lost modes, when they hit BS2. For each of these 
transformations y is transformed according to y — ► Sy, 
where S is a matrix, which is easily determined from Eqs. 
(ID), ([3]), and (|4]). The corresponding transformation of V 
is V — > SVS T . At the end of the calculation all rows and 
columns of V that represent lost modes are erased, and 
the result is a matrix, which contains the same informa- 
tion as Eqs. (fT5|) and (|16p . We have used this method to 
compute the results for the pulsed OPO in Sec. IV. 



III. FILTERING 

The OPO produces pairwise quantum correlated fields 
in a large number of cavity field modes, and in experi- 
ments it is necessary to apply a frequency filter before 
the trigger detector, see Fig. [TJ to ensure that the pho- 
ton is derived from a well defined sideband, and that the 
signal field predominantly occupies a single mode. We 
now turn to a description of such a filter modeled by a 
cavity with two beam splitters BSi and BS2 and two per- 
fectly reflecting mirrors as depicted in Fig. [3J Copying 
the notation from Sec. II, d(t) is the input field, b(t) is 
the output field, v(t) is a field in the vacuum state, Ci(t) 
is the field at different positions inside the cavity (see the 
figure), ti (rt) is the field transmission (reflection) coef- 
ficient of BS^ and r F is the round trip time in the filter 
cavity. The fields are related according to the equations 



c 1 (t) = 

h{t) = 
h{t) = 
c 4 (t) = 
b(t) = 

from which we derive 



t\d{t) + ir 1 C4(t), 
Ci(t- r F /4), 
t 2 v(t) + ir 2 c 2 (t), 
-c 3 (t-3r F /4), 
t 2 c 2 (t) +ir 2 v{t) 1 



(23) 
(24) 
(25) 
(26) 
(27) 



b(t) = ir 2 v(t) - ir 1 4^ / (r 1 r 2 ) n v(t - (ra + 1)t f ) 

71=0 

OO 

+ ht 2 ^{nr 2 ) n a{t - (n + l/4)r F ). (28) 

n=0 

Transforming the output field to frequency space 



b(uj) = ir 2 v{w) 



1 - r 1 r 2 e i " TF 



v{u>) 



1 - rir 2 e MTF 



a(w), (29) 




FIG. 3: Filter cavity. The input field a(t) is partially trans- 
mitted into the cavity through the beam splitter BSi. The 
reflected component of the input field is lost as is the part of 
the intra cavity field that leaves the cavity through BSi. The 
part of the intra cavity field that leaves the cavity through 
BS2 contributes to the output field b(t). v(t) is a field in the 
vacuum state. 



we find 



<&t( w )5( w )) = 



fZfZ 

l l l 2 



1 + r\r\ — 2r\r 2 cos(wr F ) 



(30) 

and it is apparent that the frequencies transmitted most 
efficiently through the filter are those which satisfy the 
condition tur-p = 27m, n G Z. The free spectral range of 
the filter is thus 27t/t f , and for a given r F the bandwidth 
is determined by t\ and t 2 . 

In order to select just one of the frequency modes 
emerging from the OPO, the bandwidth of the filter 
should be much smaller than the free spectral range of 
the OPO, but larger than the bandwidth of the OPO. 
Also, the free spectral range of the filter should be large 
compared to the bandwidth of the parametric down con- 
version. Experimentally it can be difficult to build cavi- 
ties that satisfy the last condition, and one may use in- 
stead a sequence of different filter cavities @. Since we 
assumed an infinite bandwidth of the down conversion in 
Sec. II, we must, for consistency, use an infinitely small 
filter cavity in the theoretical treatment. We thus assume 
(lu — uj )tf <C 1, t\ <C 1, and t 2 -C 1 for fixed K\ = tl/rp 
and k 2 = t^/r-p, where luq = 27rn/r F , n £ Z, is a reso- 
nance frequency of the filter cavity, and expand Eq. (|29|) 
to lowest order: 



b(u) =i(l- — —, , 



v(lu)+ 



a(u). (31) 



In this limit the filter transmission function is a 
Lorentzian 



(St( w )S( w )) 



(ki + k 2 ) 2 + 4(lu - lu q ) : 



(a^ (uj)d(uj)). 



(32) 



6 



We consider a single mode h(t) of the output field, and 
define h a (t) and h v (t) according to 

J h*(t)b{t)dt = J h* a (t)a(t)dt + J h* v (t)v{t)dt. (33) 

Since the vacuum state does not contribute to normally 
ordered expectation values, we shall not need h v (t) in the 
following, but from Eq. (|28|) 

h(t')^r 2 e^ oi - t '' t) -^ + ^-) (t '- t) dt l , 

(34) 

i.e., the action of the filter is effectively to transform the 
mode function h(t) into the (not properly normalized) 
mode function h a (t). Note that if the time dependent 
part of the phase of h(t) is chosen as e~ luJat , the time de- 
pendent part of the phase of h a (t) is also given as e~ luJot . 
When we choose the resonance frequency of the filter to 
equal half the carrier frequency of the OPO pump beam, 
we may thus continue to use real mode functions and omit 
the factor e luJ °^ ' ~ l) in Eq. (35). For a given bandwidth 
Ki + k 2 of the filter, 

J \K{t)\ 2 dt = 1 - J \h v (t)\ 2 dt = 

J («i + k 2 ) z + 4(lj - uj a y 

is maximal for k\ = k%, and we thus assume Ki — k 2 — K 
in the rest of the paper. 

IV. MODE FUNCTION OPTIMIZATION FOR 
HERALDED GENERATION OF A 
NON-GAUSSIAN STATE OF LIGHT 

With the necessary tools at hand we can now proceed 
to an analysis of the preparation of non-Gaussian light 
states that can be stored in atomic samples. Gaussian 
states have Wigner functions that are positive for all ar- 
guments, and as a measure of non-Gaussian character we 
shall refer to negative values of the Wigner function oc- 
curring at the origin of phase space, e.g., for odd number 
states and odd Schrodinger cat states. In Ref. [ll| we 
derived an expression for the value of the Wigner func- 
tion at the origin of the state of an arbitrary real mode 
of the multimode state generated when conditioning on 
a photo detection event 

, ]= VsMVll + V22 - 2) - V 33 V 2 \ - ViiV^ 

A ' ' Tr(V 33 V ii y/*(V 11 + V 22 -2) 

(36) 

Vjk are the elements of the Gaussian covariance matrix 
of the mode in which the APD detection takes place 
(quadrature variables 1 and 2) and the chosen mode of 
the output state (quadrature variables 3 and 4) before 
conditioning. V is computed from the definition given 
just below Eq. ([^1) by use of the mode functions of the 



two modes, the transformation (|34p of the trigger mode 
function due to the filter, Eq. (TT7| (with b(t) replaced by 
the relevant linear combination of b{t) and the annihila- 
tion operator of the vacuum field entering into the system 
at the beam splitter in Fig. [T|) , and the two-time correla- 
tion functions given in Eqs. (fl~5f and (fT6)) . We take the 
trigger mode function to be constant in a time interval 
of duration At t positioned at time U and zero otherwise, 
where Att is much shorter than all other time scales in 
the system. In Appendix A we discuss the applied de- 
tector model and the choice of trigger mode function in 
more detail. 

As explained in the Introduction the mode that will 
be stored is determined by the shape and timing of the 
strong storage pulse, and it is advantageous if one can 
initiate the generation of the storage pulse before it is 
known whether a trigger detection event will actually 
take place at the right time relative to the pulse stored. 
The protocol is thus probabilistic. We consider an at- 
tempted storage as successful if a trigger detection event 
takes place within a predefined time interval of duration 
T, whose position on the time axis is determined relative 
to the pump pulse, if the OPO is pulsed, and relative to 
the storage pulse, if the pump field is time independent. 
In practice, T is at least as large as the temporal resolu- 
tion of the APD detection system, which is of order 1 ns 
[Tot . Since the shape of the storage pulse is independent 
of the actual trigger detection time, the quantity we op- 
timize is the mean value of the Wigner function at the 
origin 

T/At t . T/At t 

W(0,0)= W<(0,0)/ J2 p i> ( 3? ) 

i=l ' 3=1 

where the sum is over the T/Att trigger modes inside the 
acceptance interval T, P, is the probability to obtain a 
detection event in trigger mode i, which for infinitesimal 
Att is equal to the expectation value of the number of 
photons in the ith trigger mode, and Wi(Q, 0) is the value 
of the Wigner function at the origin of the chosen mode of 
the generated state conditioned on a detection in trigger 
mode i given in Eq. (|36p . The total trigger probability 

T/At t 

P= Y, p i (38) 
t=i 

is small, and hence we neglect the possibility to have two 
or more detection events within the time interval T. 

Under realistic experimental conditions there will be 
significant losses, and these are included in the calcula- 
tions presented below. The chosen experimental parame- 
ters are given in Tablefl] and for the case of a pulsed OPO 
we assume that the pump field has a Gaussian envelope 
and write z(t) as 

* ) = ^r exp (-2ifJ' (39) 
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Quantity 


Symbol 


Value 


Total OPO cavity length 


L 


81 cm 


Round trip time in the OPO cavity 


r 


2.7 ns 


Loss in cavity 


T 

ri 


0.004 


Transmission of OPO output mirror 


77 1 

ti 


0.127 


Reflectance of beam splitter in Fig. ^ 


R 


0.05 


Filter bandwidth 


2k 


5.63 • 10 s s" 1 


Trigger channel efficiency 


Vt 


0.07 


Signal channel efficiency 




0.70 



TABLE I: Parameters used in the numerical minimization 
of W(0, 0). r\ t includes both propagation efficiency in the 
trigger channel and the efficiency of the APD, while rj s is 
the propagation efficiency of the signal channel. The values 
are chosen in accordance with the numbers given in Refs. 



Optimizing (|57|) for different values of T p , s, and T (or z 
and T for the case of a time independent pump field) by 
optimizing the shape of the mode function of the mode 
that is to be stored, we obtain the values of W(0, 0) 
shown in Fig. 2^, and the corresponding success prob- 
abilities are given in Fig. 0Jd. Examples of optimized 
mode functions are plotted in Figs. [5] and [Sj We note 
that the absolute temporal position of the interval T has 
been chosen to maximize the success probability. Also, 
the values of z for the OPO driven with a time indepen- 
dent pump field have been chosen in order to obtain a 
total flux of photons in the degenerate mode of the out- 
put field from the OPO, which is of order 2 ■ 10 6 s _1 as 
in Ref. [9|. This corresponds to a success probability of 
1.9 • 10~ 4 for T/t = 10. The values of s have been ad- 
justed to obtain comparable success probabilities in the 
calculations for pulsed pump fields. 

Considering the trends in Fig[4}), we find that the suc- 
cess probability increases when z increases and when s 
increases for fixed T p as expected. The success probabil- 
ity is also an increasing function of T, because more terms 
are included in the sum in Eq. (|38|) . when T increases. 
For a time independent pump field the increase is linear, 
because all the Pi 's are equal, while the success probabil- 
ity levels off to a constant value for a pulsed pump field 
when T increases beyond the temporal width of the in- 
tensity distribution of the generated output field, which 
is roughly of order the width of the intensity distribution 
of the pump pulse plus the mean lifetime of a photon 
in the OPO cavity (« 7.6 r) plus the mean lifetime of a 
photon in the filter cavity (« 0.7 r). The constant value 
is thus reached faster for short pulses. It is apparent from 
the figure that a decrease in T p must be accompanied by 
an increase in s in order to keep the success probability 
unchanged. Both of these changes increase the required 
peak value of z, i.e., a larger field strength or a larger 
conversion efficiency in the crystal is needed. 

To explain the trends in Fig.[4ji we note that W(0, 0) of 
a state with density matrix p = X)^Lo E™=o c n,m\n){m\ 
1S Yl7Lo{ c 2i,2i-C2i+\,2i+i) I and thus an optimization of 
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0.01 i — r 




5 10 15 20 25 30 



T/x 

FIG. 4: a) W(0,0) and b) success probability as a function 
of T. Solid lines: time independent continuous-wave (CW) 
pump field with z = 0.010, 0.014, 0.017, and 0.020; dashed 
lines: pulsed pump field with T p /t — 8 and s — 0.03, 0.04, 
0.05, and 0.06; dotted lines: pulsed pump field with T p /r = 3 
and s — 0.04, 0.05, and 0.06; dash-dotted line: pulsed pump 
field with T p /r = 1 and s = 0.07; dash-dotted line labeled 
SP: single pass (SP) operation for a pulsed pump field with 
Tp/r = 3 and s = 0.05. Within each series a smaller s or z 
corresponds to a more negative value of W^O, 0) and a smaller 
success probability. 



W(0,0) is equivalent to a maximization of the odd pho- 
ton number components. An ideal single mode squeezed 
vacuum state is a superposition of even photon num- 
ber states, and if a photon is annihilated from the state, 
these are converted into odd photon number states. The 
present experiment deviates from this ideal situation in 
two respects. Firstly, the generated state is a multimode 
state, since all the generated photon pairs do not belong 
to the same mode. The result is that the overlap between 
the optimal mode and the modes of the generated pho- 
ton pairs is only partial, and this introduces even photon 
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number components into the state of the optimal mode. 
The mechanism is particularly severe if photon pairs are 
generated in the outer regions of the optimal mode such 
that the overlap between the optimal mode and the mode 
of the photon pair is significantly smaller than unity, but 
also significantly larger than zero. This suggests that a 
short pump pulse, where the down conversion process is 
only turned on in a short period, will lead to a more 
negative value of W(0, 0), and this is also what is ob- 
served in Fig. rjj,. The effect is particularly pronounced 
for very large T. In this limit W(0, 0) approaches a con- 
stant value for pulsed pump fields, because Pi is small 
for trigger modes far from the center of the intensity dis- 
tribution of the generated state and these modes do thus 
not contribute significantly to the sums in Eqs. (|37|) and 
(|38p . For a time independent pump field, on the other 
hand, 1^(0,0) approaches the value of the unconditional 
state, because the success probability approaches unity 
for very large T, and all storage attempts are accepted 
as successful. The second deviation from ideal behavior 
is that losses degrade the odd photon number states into 
both odd and even photon number states, and since a 
photon number state with a large number of photons 
is more fragile than a single-photon state, it is to be 
expected that a large intensity of down converted pho- 
tons will lead to a less negative value of W(0, 0). Short 
pulses must be intense in order to keep the success prob- 
ability close to 10~ 4 , and this explains why we do not 
obtain more negative values of VF(0,0) than those pre- 
sented in Fig. 0^ if the pulse duration is decreased below 
T p /t — 1 but rather observe less negative values for very 
short pulses. 

In Fig. |31 we have included an example of single pass 
down conversion for T p /t = 3, i.e., we have increased t\ 
to unity without changing the rest of the parameters in 
Table |TJ When the OPO cavity is absent, the temporal 
width of the light pulse is only broadened in the filter 
cavity, and this broadening is small since the mean life- 
time of a photon in this cavity is only 0.7 r. It is thus 
possible to determine from the precise APD photo detec- 
tion time whether the photon pair that gives rise to the 
APD detection was generated in the beginning, in the 
middle, or in the end of the pulse. For small T this leads 
to a negative value of W(0, 0) because we know rather 
precisely where the second photon in the pair is, if it has 
not been lost. On the other hand, if we do not discrim- 
inate between early and late trigger detection events a 
single optimal compromise will only have a small overlap 
with the actual modes, which are well localized in time. 
The resulting rapid increase towards positive values in 
W(0,0) with increasing T can be avoided by increasing 
the mean lifetime of a photon in the filter cavity, corre- 
sponding to a smaller frequency width of the filter, and 
thus the intensity of the light transmitted through the 
filter will decrease, which leads to an even smaller suc- 
cess probability than in Fig. [Dd. The success probability 
can be increased by increasing s, but this will also result 
in a less negative value of W(0, 0), and from these argu- 



a) 

0.35 




05 1 1 1 1 -0.05 1 ' 1 1 

-30 10 50 90 -30 10 50 90 

tlx 111 



FIG. 5: Mode functions for the modes of the generated state 
with the most negative values of W(0, 0) for a pulsed pump 
field with Tp/r = 8 and s = 0.05. t = is chosen to coincide 
with the time where z is maximal, a) T/t — 0.4 (the inset is 
an enlarged view), b) T/t = 2.4, and c) T/t = 10. 



ments it is preferable to apply pulsed fields and the OPO 
cavity in the setup to get the most efficient generation of 
non-Gaussian states for storage. 

The values of W(0, 0) presented in Fig. [4ji are signifi- 
cantly above the theoretical minimum of and we 
note that the results are all obtained for physical param- 
eters, which are already achieved in experiments. More 
negative values of W(0, 0) can be reached if it is possi- 
ble to reduce the losses. If, for instance, we increase the 
signal channel transmission from 0.7 to unity, W(0, 0) is 
decreased to —0.23 for T p /t = 3, s = 0.05, and large 
T, while the success probability is unchanged, and if we 
further assume zero loss in the cavity and reduce R to 

0. 01, we find W(0, 0) = —0.27. An increase in the trigger 
channel transmission or in the APD detector efficiency 
will increase the success probability, and this increase can 
be transformed into a more negative value of W(0, 0) by 
decreasing z, s, T, or R. 

Considering the optimized mode functions for a pulsed 
pump field in Fig. O it is apparent that a fast variation 
in the mode function is present for small values of T, but 
this variation is smoothed out for larger T. The inset in 
Fig. [5^, shows that the period of the fast variation is r, 

1. e., the round trip time in the OPO cavity. A photon 
generated in one of the time localized modes in the OPO 
cavity can only leave the cavity at times separated by 
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FIG. 6: Pump pulse (dash-dotted line) and mode function 
of the mode of the generated state with the smallest value 
of W(0,0) (solid line) for T p /r = 3, s = 0.05, and T/r = 
32.4. The dotted line is the function in Eq. (|40)) , while the 
approximation (|41[) (dashed line) is almost indistinguishable 
from the optimum mode (solid line) in the figure. W(0, 0) = 
-0.0591 for the dotted line, W(0, 0) = -0.0610 for the dashed 
line, and W(0,0) = -0.0611 for the solid line. 



an integer number of r's, and thus, in the output field 
from the OPO, the two photons in a pair are separated 
by an integer number of r's. The fast variation in the 
optimal mode function can thus be explained from the 
fact that the mean lifetime of a photon in the filter cav- 
ity is smaller than r, since this means that the observed 
photon in a pair is typically delayed less than r in the 
filter cavity, and the detection time provides partial in- 
formation on the time of generation of the photon pair. 
The information is erased as T approaches r, and this 
leads to the steep rise in W(0, 0) towards less negative 
values observed for all curves to the very left in Fig. [4ji. 

If the duration of the pump pulse is short compared to 
r, the optimal mode function differs qualitatively from 
those presented in Fig. [5] In this case the optimal mode 
function is a series of spikes separated by r, and the width 
of the spikes is determined by the width of the pump 
pulse provided the width of the pump pulse is also short 
compared to T. In the opposite limit of a time indepen- 
dent pump field, the optimal mode function for short T 
is qualitatively a function, which has a maximum close to 
the time of the trigger detection and decays exponentially 
when moving to the right or left from this maximum. Su- 
perimposed on this is a fast variation with period r. The 
mode function is not completely symmetric around the 
time of the trigger detection due to the filtering. For 
larger values of T the fast variation disappears, and the 
maximum becomes more rounded. 

In Fig. [S] we show the temporal shape of the pump 
field and the optimal mode function for T/r = 32.4 in 
the same graph, and it is apparent that the generated 



output field is both broadened and delayed due to the 
nonzero lifetime of a photon in the OPO cavity as ex- 
pected. Since W(0, 0) and P are both independent of T 
for large T, and since we have to choose the same mode 
function independent of the precise time for the APD de- 
tection, the optimal mode function is determined solely 
from the generated output field when T is large, and it 
seems reasonable that the optimal mode function just fol- 
lows the square root of the intensity distribution of the 
generated state. Without applying the above theory, a 
crude guess for the optimal mode function is 



h{t) 



z(Oexp^-^±i(t-t'))^', 



(40) 



and this function, which is also shown in the figure, is 
actually quite close to optimal. With the above theory at 
hand, however, it is easy to calculate the actual intensity 
distribution from Eq. (|15p . and this leads to the mode 
function 



h(t) (x 



\ 



fn+l 



^(rii 2 ) 2 "sinh 2 ^ z{t - kr) 



(41) 



\k=l 



which is even closer to the optimal mode function. Due 
to its simplicity, Eq. (|4"Tj) can be useful in a storage ex- 
periment, where T must be sufficiently large to obtain a 
satisfactory success probability. We note, however, that 
neither Eq. (|4H|) nor Eq. (|4Tj) is able to reproduce the 
optimal mode function in Fig. [5^ for T/r = 0.4. 



V. CONCLUSION 

In conclusion, we have presented a multimode descrip- 
tion in time domain of the optical parametric oscilla- 
tor, which is valid for both pulsed and continuous- wave 
pump fields, and we have used the theory to analyze 
non-Gaussian states that can be stored in atomic sam- 
ples. So far light storage has been demonstrated with 
Gaussian states only, and the storage and retrieval of a 
non-Gaussian state is a hall mark in the demonstration 
of atom-light quantum interfaces. Our analysis suggests 
that, with similar success probabilities, pulsed pump 
fields lead to more negative values of the Wigner function 
of the stored state than time independent continuous- 
wave pump fields. 

An essential ingredient in the experiment is the tem- 
poral spreading induced by the OPO and filter cavities, 
because it separates the two photons in a generated pho- 
ton pair temporally by an unknown amount and ensure 
that it is impossible to infer the precise position of the 
second photon in a pair from the time at which a detec- 
tion took place in the APD. Without this spreading we 
do not obtain large negative values of W(0,0) for large 
T. After having thus characterized the system and iden- 
tified the optimum strategy for given setup and pulse 
parameters, one may now move a step further and try to 
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design filters and pump pulses in order to generate opti- 
mal mode functions which have high non-Gaussian state 
content and which are particularly easy to handle in the 
storage part of an experiment. 

In the present treatment we have chosen to use the neg- 
ative value attained by the Wigner function as a measure 
of the non-classical character of the state generated, but 
we note that only minor changes are required in order to 
optimize other features as for example the single-photon 
or Schrodinger kitten state fidelity. 

The authors acknowledge discussions with Jonas S. 
Neergaard- Nielsen, Anders S. S0rensen, and Eugene S. 
Polzik. 



APPENDIX A: DISCUSSION OF THE APPLIED 
DETECTOR MODEL 

Various detector models are discussed in the literature. 
Considering first detection in a single mode, a detector 
may for instance be a photon number resolving detector 
or an on/off detector. The photon number resolving de- 
tector registers the number of photons in the detected 
mode, and the detected mode is projected onto the de- 
tected photon number state, while the on/off detector 
only distinguishes between the outcomes 'vacuum' and 
'not vacuum', and, depending on the measurement re- 
sult, the detected mode is projected onto one of these 
subspaces. If a continuous beam of light is observed, a 
third detector model is often assumed, where a detection 
event at time t corresponds to an annihilation of a pho- 
ton in a mode of infinitesimal temporal width positioned 
at time t. Eq. <|36j) is derived under the assumption that 
a detection event is equivalent to application of the anni- 
hilation operator of the observed mode to the state of the 
system, but since the flux of trigger photons is small in 
the experiment considered in Sec. IV and the possibility 
to have two or more photons in a trigger mode can be 
neglected, all three detector models lead to practically 
identical results. 

A real detector consists of a medium, which absorbs 
photons, and typically several different microscopic tran- 
sitions in the medium are allowed. The detector is then 
able to detect photons in more than a single mode. In an 



APD a photon absorption is accompanied by an excita- 
tion of an electron from a bound state to a free or a solid 
state conduction band state, the signal is amplified, and 
the resulting photo current is measured. The informa- 
tion concerning the precise microscopic transition is lost 
in the amplification process, and, when we condition on 
the macroscopic detection event, we must average over 
all modes that could have led to the observed detection 
(see Since the filter in front of the detector se- 

lects a narrow band of frequencies, we assume that the 
detector is equally sensitive to all frequencies transmit- 
ted by the filter, i.e., the detector bandwidth is assumed 
to be infinitely broad compared to the filter bandwidth, 
and all modes are observed by the detector. In this case 
we should sum over a complete set of modes on the in- 
terval T in Eq. (JSTJ) . The obtained values of W(0,0) 
and P are independent of the choice of basis which may 
be seen from the following argument. Assume that the 
mode functions f+ and the mode functions gi both con- 
stitute a complete basis on the time interval T. fa and gi 
are related by a unitary transformation = y\ XJ^g^ 
which implies that the corresponding mode annihilation 
operators a{ = j f*(t)a(t)dt and af = j g*(t)a(t)dt are 
related by aj = £V JJ*jO,j. If a photon is annihilated in 
the mode a{, the density matrix is transformed accord- 
ing to p — > a{ p(a{ )t/Tr(a{ p[a{ where Tr denotes the 
trace. To perform the averaging we multiply by the prob- 
ability Tr(af p(a{ )t) to annihilate a photon in the mode 
aj , and sum over all i. After normalization this leads to 



Tr(yVafp(a{)t) 



(Al) 



Since Y,i°>iP(H) 



N = 



J2j J2k ®j Pi^kV^jk = J2i 6? p(af ) f , the transformed den- 
sity matrix is independent of the choice of basis, and the 
Wigner function is then also independent of the choice of 
basis. It is thus not necessary to use the modes, in which 
the detections actually take place, and for mathematical 
convenience, we have chosen to sum over time localized 
modes in Eq. (|37|) . 
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